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I.  INTRODUCTION 

The  purpose  of  this  technical  memorandum  is  to  describe 
an  analytical  method  for  obtaining  the  sound  pressure  field  due 
to  a  plane  wave  incident  on  a  complex  consisting  of  an  infinite 
length  elliptical  ("nearly  straight")  baffle  and  an  infinite 
length  right  circular  cylinder.  The  problem  is  formulated  so 
that  the  baffle  shields  the  cylinder  from  the  incident  plane 
wave.  The  model  used  for  this  boundary  value  problem  is  a  gen¬ 
eralization  of  a  previous  concentric  baffle  model  used  for  studying 
the  interaction  of  the  sonar  baffle  and  transducer  with  screw 

i  p 

noise.  ’  The  analytical  results  are  presented  in  a  form  suitable 
for  computing  the  total  sound  pressure  at  the  transducer 
(cylinder)  face  as  well  as  the  incident  and  reflected  portions 
of  the  field  solution  there. 

The  material  in  this  memorandum  is  devoted  solely  to 
the  development  of  a  formal  solution  for  the  boundary  value 
problem  as  described.  An  investigation  of  numerical  methods  for 
the  computation  of  the  solution  is  being  conducted.  In  addition 
it  should  be  noted  that  while  the  boundary  value  problem  described 
in  this  memorandum  is  useful  for  studying  the  interaction  of  the 
baff le-transducer  with  screw  noise,  the  mathematical  methods 
which  are  developed  can  readily  be  applied  to  a  study  of  the 
interaction  of  the  transducer  and  an  elliptical  dome  during 
transmission.  This  application  is  currently  being  studied. 


n 

'‘‘An  Analytic  Solution  to  the  Sound  Pressure  Field  Resulting  from 
a  Plane  Wave  Incident  on  a  Cylinder  and  Concentric  Cylindrical 
Section,  TRACOR  Report  63-II2-U,  March  14,  1963 . 

O 

An  Analytical  Investigation  of  Sonar  Dome  Spoking,  TRACOR 
Report,  64-llo-C,  February  14,  1964. 
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II.  FORMULATION 

In  the  formulation  of  the  solution  of  our  problem  we 
consider  the  four  coordinate  systems  shown  in  Figure  1.  The 
sound  pressure  as  a  function  of  the  space  coordinates  and  time 
satisfies  the  wave  equation 


* 

xx 


+ 


(1) 


The  subscripts  in  this  equation  denote  partial  differenta tion, 
the  partial  derivatives  being  taken  with  respect  to  the  entity 
used  as  the  subscript.  The  equation  may  be  separated  by  assuming 
that 


f(x:y,t)  =  p(x,y)  f(t)  .  (2) 


Substitution  of  this  assumed  form  for  'I '  into  equation  (1)  yields 


where  the  dot  notation  represents  time  derivatives.  The  left 
hand  side  of  equation  (3)  is  a  function  of  x  and  y  while  the 
right  hand  side  is  a  function  of  time  only.  Since  x,  y,  and 
t  are  independent  each  side  of  equation  is  a  constant.  Thus 
taking  the  constant  (which  is  called  the  separation  constant) 

p 

to  be  -k  equation  (3)  yields 


(4) 

(5) 

(6) 
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This  is  the  Helmholtz  equation  and  it  is  the  partial  differential 
which  describes  the  spatial  variation  of  the  sound  pressure  field. 

■  ”  The  coordinate  systems  shown  in  Figure  1  are  related  by 


y  =  l  +  x'  j 

x  -  -y'  ,  (7) 

x  =  75  cosh  [i  cos  cp 

x  =  fj  sinh  \x  sin  cp  ,  (8) 

and 

x !  =  p  cos  0 

y'  =  p  sin  G  .  (9) 


When  the  Helmholtz  equation  is  transformed  into  elliptical 
coordinates  through  the  equations  (8)  it  takes  on  the  form 


2  2 

pM-|i  +  pcpcp  +  H"  £cosh2  M-  -  cos2  cp]p  =  0  .  (10) 

This  form  of  the  Helmholtz  equation  can  be  separated  by  assuming 
a  solution  of  the  form 


p  =  M(n)  P(cp)  . 

The  resulting  ordinary  differential  equations  are 

Pn  +  (b  -  h2  cos2  cp)P  =  0 

and 

M,!  +  (h2  cosh2  p-  -  b)M  =  0  . 
Here  b  is  the  separation  constant  and 


(ID 


(12) 

(13) 


(14) 


The  equations  (12)  and  (13)  are  forms  of  the  Mathieu  equation. 
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At  this  point  a  few  remarks  regarding  solutions  to  the 
Mathieu  equations  are  appropriate.  We  are  attempting  to  define 
a  sound  field  continuous  near  the  baffle  (|JL=e)  .  This  implies 
that  we  seek  solutions  to  equation  (12)  which  are  periodic  in  cp 
(where  2n  is  an  integral  multiple  of  the  period) .  There  exist 
two  countably  infinite  sets  of  values  for  the  separation  constant 

.  .,  and  (3q,  P2  .  ..,  for  which  the  solutions 
.  In  the  problem  at  hand  we  are  interested 
only  in  those  solutions  which  are  even  functions  of  cp  about 
cp  =  tt/2.  In  this  particular  case  we  need  only  study  the  solutions 
which  correspond  to  the  subsets  <Xq,  a 2,  . ..,  and  p^,  p^,  P5  •••> 

since  these  and  only  these  have  this  desired  property,  i.e. 
are  even  functions  of  cp  about  cp  =  tt/2. 

The  above  requirement  for  a  specific  set  of  solutions 
for  equation  (12)  defines  the  allowable  values  of  b.  For  each 
value  of  b  there  are  two  solutions  of  equation  (13) • 


Solutions  of  equation  (12)  correspondong  to  the  a q,  a-^, 
a2  . . .  are  denoted  by  Sem(h,cos  cp)  and  those  corresponding  to 
Pq,  Pg  •  •  «  are  denoted  by  Som(h,cos  cp)  . 

The  two  solutions  of  equation  (13)  corresponding  to  the 
am  are  denoted  by  Jein(h,cosh  p.)  and  Nem(h,cosh  p.)  ,  whereas  the 
two  corresponding  to  the  Pm  are  denoted  by  Jom(h,cosh  p.)  and 
No  fh, cosh  p.) . 

We  therefore  have  the  solution  for  the  F^lmholtz  equation 
which  is  an  even  function  of  cp  about  cp  =  tt/2, 


P  = 


'  2m  Je2m  +  B  2m  Se2m 


*  ^ '  2m-rl  Jo2m+l  +  D  2m+l  No2m+l 


]  So 


2m+l 


(15) 


^Morse,  P.M.  and  K.  Feshbach;  Methods  of  Theoretical  Physics; 
McGraw  Hill  Book  Company,  Inc.,  New  York,  1953;  p.  554. 
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The  expression  for  the  pressure  due  to  a  plane  wave 
proceeding  in  the  direction  of  the  positive  y  axis  (with  time 

«,  f*  i|- 

factor  e  surpressed)  in  terms  of  the  Mathieu  functions  is 


=  2  I  (■ 

m=0  L 


Se-  (h,0) 

1)  -  Se~(h,cos  cp)Jep^(h,cosh  p) 

Me0m(h)  **  ** 


+  i(-l) 


m  So2m+Zh’°) 


^°2m+l (h>cos  <?>)^°2nri*l Oi>eosh  M-)  >  •  (l6) 


Here  the  and  M°2nrt-1  (h)  are  constants. 

The  total  sound  pressure  field  is  the  sum  of  the  plane 
wave  pressure  and  the  pressure  given  in  equation  (15),  i.e. 


/8tt 


_y  r 

TT  Zi  S 


m  SeOm(h,0) 

7TT>Je2ni^h,cosh  *> 
M2m+l(h) 


+  Ne^^jcosh  n)3  Se^Ch,  cos  cp) 


iC-D^Oj.,  (h,0) 

+  f  (c2mfi  +  — x; — jr: - ^(h.cosh^ 

M2ort-l(h) 


+  D2nrt-lNo2m+l(h>cosh  ^t)^So2nrt'l (h>cos  <P>  j  (17) 

where  /fa  A ^  =  A^_,  /fa  =  Bgm,  /BnCgm+l  “  c  ^nri-i  an<i 

^  D2m+1  =  D2m  - 


hIbid,  p.  1410. 
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We  impose  the  condition  corresponding  to  zero  pressure  at 
the  baffle 

ptot^e,cp)  =  0  (l8) 

on  the  total  pressure.  This  shows 


(-l)mSe2m(h,0)^  Je^ChjCosh  e) 

+  ^  NeTThTcosFTy 


B2m  “  “  (A 


2m 


Me2m(h) 


'2m' 


(19) 


and 


i(-l)m  SOo., (h, 0)  (h,cosh  e) 

+  - - — -  .  (20  ) 


D2m+1  “  ‘  (C2m+l 


>Wh> 


No 


2ro+l (hjCosh  e) 


Equation  (1”)  is  therefore  reduced  to 


_  V  J  (“1)  SePm(h>0) 

ptot  -  ^  L  \  fA2m  +  - - HJe^h.cosh  p) 


m=0 


Je2m(hJcosh  e) 

'  Ne2m(h,  cosh  e)  Ne2m(h<  cosh  M-)]Se2m(h>  cos  <P> 


+  [C2m+l  + 


tr 


2rrH~l 


(h) 


]CJo2m+l(h»  cosh  V) 


Jo2m+i(h>  cosh  e) 
N°2n>fi  (h>  c°sh  e) 


No^^iCh,  cosh  MOlSOg^Ch,  cos  <p)l  .  (21) 


The  total  pressure  as  given  by  the  equation  (21)  represents  the 
collection  of  all  fields  which  are  even  functions  of  x,  i.e.  about 
the  y  axis,  with  zero  pressure  at  the  baffle,  and  in  which  a  plane 
wave  is  moving  in  the  positive  y  direction.  If  we  could  impose  the 
boundary  condition  at  the  face  of  the  transducer  we 
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could  complete  our  problem.  This  however  is  not  plausible  because 
of  the  nature  of  the  transducer  face  as  described  in  the  |i,cp 
coordinates . 

The  Helmholtz  equation  when  transformed  into  the  polar 
coordinates  is 


p  +  - 

PP  P 

PP 

f  h  Pee  +  k2P  =  °- 

P 

(22) 

Separation  of 

equation  (22) 

by 

taking 

P 

=  R(p)  T(0) 

(23) 

gives 

O  2 

R"  + 

-  R 
o 

’  +  (k2  -  ~)R  =  0 

(24) 

and 

r' 

P 

ipl  * 

+  n2  T  =  0  . 

(25) 

Continuity  in  0  implies  that  n  must  be  an  integer.  The  solution 
of  equation  (24)  is 


R  -  En  H^Ocp)  +  Fn  H^>  (kp)  ,  (26) 

where  the  (kp)  and  (kp)  are  Hankel  functions.  Equation 

(25)  is  the  harmonic  equation  with  sin  n0  and  cos  n0  as  solutions. 

The  solution  of  equation  (22),  which  is  an  even  ^unction 
of  0  about  0  =  0  is 


=  Z  CEn  HnX)  (kP>  +  Fn  H,$2)  (kP>3  cos  n0  •  (27) 


n=0 


THp  2  CjirO  HnO  f* n  r\1  ono  f.sp tro  £  o r*  /l  ^  \  -?  M 

—  w  v.v*'-  ww  wnv.  vVu  v  v-  U x.  cs^  «  y  Xtt  pUiai 

coordinates  is 
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CO 


pw 


2  elU  (£n^n)  (^)  (Up)  +  H^2)  (kp)]cos  ne  (28) 


n=0 


where  e  =  i  n  2 
n  2  n  >0 


The  total  pressure  is 


ptot  -  2  {<En  +  eikt  <kP) 

n=0 


ikt  en1 


.n 


r(2) 


1 


+  (Fn  +  e  ^  ~ — )H^;  (kp)  )  cos  nB 


(29) 


If  we  impose  the  condition  that  the  particle  velocity  on  the 
transducer  face  is  zero  we  find 

lv,  ein  ,  €  in  '  (kr) 

F  +  elU  -fa-  =  -  (E  +  3iU  -V)  -fan - —  . 

2  (kr  ) 


(30) 


The  expression  for  total  pressure  is  therefore  reduced  co 


eninKH^1)(kp) 


HnX) ' (kro> 
^2)  '  <kro> 


H<2)(kp)] 


cos  n0 . 


(31) 


Equation  (31)  represents  the  collection  of  all  fields 
even  about  the  y  axis  with  zero  particle  velocity  at  the 
transducer  and  in  which  a  plane  wave  is  moving  in  the  positive 
y  direction.  As  before,  if  we  could  impose  the  boundary  con¬ 
dition  at  the  baffle  on  equation  (31)  our  problem  would  be 
complete.  This  is  not  plausible  because  of  the  nature  of  the 
equation  for  the  baffle  face  in  polar  coordinates . 
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Equations  (21)  and  (31)  are  'almost"  complete 
representations  for  the  sound  field  for  the  on  axis  plane  wave 
incident  on  a  soft  elliptical  baffle-hard  transducer  complex. 
Each  of  the  collections  of  functions  contains  a  valid  representa 
tion  of  the  solution  function  in  some  domain.  Fortunately  these 
domains  overlap.  This  follows  from  the  fact  that  the  solution 
function  is  analytic  in  the  free  finite  field.  It  is  therefore 
possible  to  determine  the  and  in  equation  (21)  and  Er 

in  equation  (31)  by  "matching"  the  two  series  at  a  point  between 
the  baffle  and  transducer. 
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III. 


MATCHING  THE  TWO  REPRESENTATIONS 


In  order  to  match  the  two  representations  we  choose 
arbitrarily  the  point  cp  =  tt/2,  p.  =  (3  with  e  <  {3  <  arc  sinh  [—  ], 

or  in  polar  coordinates  0=n,  p^=[-t-a/2  sinh  [3]  for  the 
matching.  At  this  point  we  will  rnat^h  the  function  and  all 
partial  derivatives  (note  that  all  odd  derivatives  in  the  x 
direction  are  zero)  to  get 


2q+s 

5  M  Pj- 


tot 


[p(x',y'),  0(x',y')] 


a2q+s  P 


tot 


[>(x,y},  cp(x,y)] 


dy 


2q  3X'S 


ax2q  ays 


(32) 


at  0  =  tt,  p  =  p-^,  M-  =  (3,  and  cp  =  n/2 
for  q  =  0,1,  2  ...  and  s  =  0 ,  1 ,  2  . . . 


The  equations  (32)  form  an  infinite  set  of  non- 
homogeneous  linear  equations  in  the  as  yet  undetermined  constants 
^2m5  E2mJ-'l  anc*  En*  T^e  so^ution  °f  this  system  and  hence  the 
description  of  the  sound  field  may  be  obtained  to  any  desired 
accuracy  by  truncating  the  series  expressions  for  Ptot[M'>cp] 
and  Ptot[p  j9]  at  appropriate  values  of  m  and  n  and  imposing 
on  these  constants  the  proper  finite  set  of  conditions  from 
the  set  (32)  . 
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IV.  CONCLUS  ION 

The  sound  field  near  the  transducer  is  given  by  the 
truncation  of  the  solution  to  the  Helmholtz  equation  as  given 
in  equation  (31)-  The  scattered  and  incident  portions  of  the 
field  are  given  by  those  terms  of  equation  (31)  which  involve 
Hankel  functions  of  the  first  and  second  kinds  respectively.^ 
Hence,  the  total  sound  field  and  the  scattered  and  incident 
portions  thereof  can  be  computed  using  the  technique  described 
above  . 

The  technique  used  in  solving  this  elliptical  baffle- 
circular  transducer  problem  can  also  be  applied  to  other  sonar 
problems.  For  example  the  problem  of  determining  the  interaction 
of  a  circular  transducer  inside  an  elliptical  dome  during 
transmission  can  be  treated  with  precisely  the  same  formal 
rnathematical  tools  as  are  used  in  this  problem. 


5Ibid.,  p.  1371 
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